High- Visibility Multi-Photon Interference of Classical Light 
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It is shown that the visibihty of multi-photon interference for classical sources grows rapidly with 
the order of interference. For three-photon and four-photon interference of two coherent sources, the 
visibility can be as high as 81.9% and 94.4%, respectively, - much higher than the 'classical limit' 
of two-photon interference (50%). High- visibility three-photon and four-photon interference has been 
observed in experiment, for both coherent and pseudo-thermal light. 
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It is known that the nonclassical phenomena of two- 
photon interference [l| and two- photon ghost diffraction 
and imaging have classical counterparts. Two- 

photon interference of classical light has been first discov- 
ered in the pioneering experiments by Hanbury Brown 
and Twiss '4\ and since then was observed with vari- 
ous sources, including pseudothermal ones [5], true ther- 
mal ones [6], and coherent ones [7[. Somewhat later, 
ghost imaging with classical light has been demonstrated, 
both in the near-field and far-field domains [lo| . 
The only disadvantage of classical light with respect to 
two-photon interference and two-photon ghost imaging, 
compared to two-photon entangled sources, is the lim- 
ited visibility, which is always below 50% 12 1, [iS^. On 
the other hand, in interference and imaging experiments 
with classical sources the visibility is independent of the 
intensity, which can therefore be arbitrarily high. En- 
tangled sources, in contrast, should be sufficiently weak 
to provide high-visibility interference: while the visibil- 
ity is close to 100% for faint two-photon light, it drops 
with the increase in the mean number of photons per 
mode [3], 

In this paper we show that classical sources can pro- 
vide much better visibility if one passes to higher-order 
intensity correlations. This fact, which has considerable 
importance for multi-photon imaging and multi-photon 
lithography, is also remarkable in view of the numer- 
ous recent experiments on higher-order correlations and 
multi-photon interference [l5]. Indeed, while in the case 
of two-photon interference any experiment with the visi- 
bility exceeding 50% can be interpreted as non-classical, 
no such criterion is formulated for higher-order interfer- 
ence. Our results show that for multi-photon interference 
to be recognized as nonclassical, its visibility should ex- 
ceed really high values: 81.9% in the three-photon case 
and 94.4% in the four-photon one. These values are clas- 
sical visibility limits for three- and four-photon interfer- 
ence, respectively. 

For our consideration, we chose the scheme of the two- 
slit Young's interference experiment [l3|. This geome- 
try was used in many experiments on two-photon inter- 
ference and two- photon ghost imaging 0, [H], 0, 



[8|, [lfl|. This time, however, we consider the interference 
to be registered by three detectors instead of two (Fig.l), 
each detector measuring the instantaneous intensity and 
the triple photocount coincidences being counted by a 
coincidence circuit. This is the standard experimental 
technique to measure the third-order Glauber's intensity 
correlation function (ICF) lj|. Let the two sources A, B 
be classical ones, having the same statistics, the same av- 
erage intensities, and independently fluctuating phases. 
The last condition provides the erasure of the first-order 
interference in the far-field zone. However, a simple cal- 
culation shows that the intensity correlation functions 
will be sensitive to the positions of the detectors. Pass- 
ing to the normalized third-order correlation function for 
the detectors placed at points 1, 2, 3, 
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Ij being the instantaneous intensity at point j = 1,2,3, 
we obtain that 
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where f/^^-* , g*-'^^ are, respectively, the second-order and 
third-order normalized ICFs for each of the two sources 
and (t>.ij = (t>Ai - 4>Aj - 4>Bi + (t>Bj, (t>sd being the phase 
gained by the radiation of source S = A, B on the way 
to detector d = 1, 2, 3. 

According to Eq.((2l), the maximal visibility is pro- 
vided by the radiation sources with the minimal ratio 
of g'-'^Vs^^^- Among the classical sources, this 'visibility 
limit' is achieved for coherent light and is equal to 81.9%. 
Thermal radiation gives a lower visibility, 60%, which is 
still much higher than the corresponding value in the case 
of two-photon interference (33%). 

In our experiment, the third-order Glauber's corre- 
lation function was measured through the coincidence 
counting rate of three detectors (Fig. 2). As the radiation 
source, we used a frequency doubled Q-switched YAG:Nd 
laser with the wavelength 532 nm, pulse duration 5 ns, 
and the repetition rate 3 kHz. Instead of two slits, which 
should be very precisely matched in width to achieve the 
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maximum visibility of multi-photon interference, a single 
slit of width 150/i was used, followed by a birefringent 
crystal (calcite). The crystal split the beam into the or- 
dinary one and the extraordinary one; with the slit and 
the crystal placed between crossed polarization (Glan) 
prisms, this configuration was equivalent to two identical 
slits separated by a distance of 1.3 mm. In the far-field 
zone, where the interference pattern was formed, the ra- 
diation was attenuated using neutral-density filters and 
fed into a three-arm Hanbury Brown-Twiss interferom- 
eter with three photon-counting avalanche diodes and a 
triple-coincidence circuit. Attenuation was necessary to 
keep the average number of photocounts per pulse much 
less than one; otherwise, because of the detectors dead- 
time effect, the photocount and coincidence rates would 
be measured incorrectly. Gating of the registration elec- 
tronic system provided suppression of the dark noise by 
several orders of magnitude. In order to scan the in- 
terference pattern, thick glass plates were placed at the 
inputs of detectors 1 and 3. Spatial mode selection was 
provided by 70^ apertures Al, A2, A3 placed in front 
of each detector. By turning the glass plates, one could 
scan the phase of either one or two detectors within the 
range . . . 6tt. The phase of detector 2 could be adjusted 
by moving the aperture A2. 

The varying relative phase of the two sources was intro- 
duced by inserting an electro-optical modulator (EOM) 
after the calcite crystal. The EOM was oriented in such 
a way that its bias voltage induced a phase shift be- 
tween the ordinary and extraordinary beams. Apply- 
ing an AC voltage with the amplitude 85 V (slightly be- 
low the quarter- wave one) and the frequency 50 Hz, we 
erased the interference pattern in the time-averaged in- 
tensity distribution. Although this harmonic phase mod- 
ulation was different from a randomly varying relative 
phase of sources A, B (Fig.l), it resulted in the same 
time-averaged ICFs. 

We studied third-order interference both for coherent 
and pseudo-thermal radiation at the input. Pseudo- 
thermal light was prepared by means of a rotating 
ground-glass disc placed after the calcite crystal. The 
envelope of the third-order spatial ICF was determined 
by the coherence radius of the radiation, which, in its 
turn, depended on the sizes of the spots formed on the 
disc by the ordinary and extraordinary beams, and its 
FWHM corresponded to approximately two periods of 
the interference pattern. Because of this, the visibility in 
the case of thermal light was considerably less than the 
theoretical value. 

The results for the coherent and pseudo-thermal cases 
are presented in Fig. 3 a,b, respectively. Experimen- 
tal points correspond to the normalized third-order 
Glauber's IGF; curves show the dependencies given by 
Eq.(l2|). In Fig. 3b, the finite value of the coherence ra- 
dius was also taken into account. From Eq.Q one can 
see that the maximum visibility is achieved when two of 



the three phases are varied simultaneously in the oppo- 
site directions, 012 = —(f>32- In accordance with this, in 
our experiment the glass plates in front of detectors 1 
and 3 were rotated synchronously, both clockwise (since 
detector 1 was in the reflected beam and detector 3 in 
the transmitted beam, this led to the opposite variation 
of the phases). The obtained visibility for the case of co- 
herent radiation is 74%; for the case of pseudo-thermal 
radiation, 38%. In each plot, we also show the spatial de- 
pendence of single counts for one of the detectors whose 
phase was scanned. Although the dependence is not com- 
pletely fiat (the variation is caused by the speckle struc- 
ture of the laser light and the envelope of the single-slit 
diffraction pattern), the two-slit interference pattern in 
the intensity distribution is erased. Note that the pre- 
sented ICF was normalized to the product of the three 
intensities; as a result, the 'noisy' structure of single- 
photon counts did not infiuencc the third-order interfer- 
ence pattern. 

If the third-order interference pattern is scanned by 
only one of the three detectors, the visibility is lower 
than 81.9% but still considerably higher than in the case 
of the second-order interference. It reaches its maximal 
value (approximately 70.5%) when the relative phase 
between the remaining two detectors is equal to 7r/2. 

Quite similarly, one can show that the fourth-order 
intensity interference of two classical sources is given by 
the formula 
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Here, the same notation for the phases and normalized 
ICFs is used as in Eq.([2|). Analysis of this expression 
shows that the maximal visibility of the fourth-order in- 
terference for thermal sources is 77.8% while for coherent 
sources, it is 94.4%. The last figure exceeds the visibility 
values achieved in nearly all known 'four-photon' exper- 
iments. 

Experimental observation of the fourth-order interfer- 
ence by registering four-fold coincidences was difficult be- 
cause of the low coincidence counting rate in this case. 
For this reason, we turned to another method of mea- 
suring spatial ICFs, the one based on digital image pro- 
cessing (see, for instance, Ref. [13] )• The interference 
pattern in the far-field zone was registered by a digi- 
tal photographic camera Canon PowerShot S2 IS. The 
source, again a frequency-doubled YAG:Nd laser, in this 
case had a repetition rate 50 Hz. Each frame was made 
with a single laser pulse, the exposure time (1/60 s) be- 
ing less than the distance between the pulses. A typical 
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interference pattern recorded in one frame is shown in 
Fig. 4. To accumulate sufficient statistics, n = 500 shots 
were made. Due to the phase shift introduced by the 
EOM, the phase of the interference pattern varied from 
frame to frame, so that the intensity spatial distribution 
averaged over all frames revealed almost no interference 
(the visibility of Young's interference pattern was less 
than 10%). In each frame, a rectangular area was se- 
lected, with the dimensions 50 pixels in the vertical (y) 
direction and 600 pixels in the horizontal (x) direction (a 
rectangle shown in Fig. 4). First, the intensity distribu- 
tion recorded in each frame was averaged over the y side 
of the rectangle, so that intensity distributions Ij{x) for 
all n pulses were obtained. Next, averaging over j was 
performed, the averaged intensity distribution and the 
third- and fourth-order normalized correlation functions 
being calculated as 
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Note that the arguments of intensities used in (l5|), (l6|) 
are chosen so as to provide the maximal visibility. In 
the above-described coincidence method of g^^'^ measure- 
ment, this was achieved by scanning two detectors in op- 
posite directions and the third one being fixed; for a g^^^ 
measurement, the fourth detector should be scanned with 
a double speed. 

Fig. 5 shows the third-order (a) and fourth-order (b) 
interference patterns obtained for coherent sources by 
means of digital-image processing. For convenience and 
similarity with Fig. 3, the x coordinate, originally mea- 
sured in pixels, is plotted in phase units, so that x = 0i2; 
the other phases obey the relations = —4>i2,4>A2 = 
— 2(/)i2. As expected, the distribution in Fig 5a is simi- 
lar to the third-order interference pattern registered by 
means of coincidence method. The interference visibili- 
ties achieved for three- and four-photon interference are 
73% and 93%, respectively. For comparison, theoretical 
dependencies given by Eqs ^ , ^ are shown as thin solid 
lines. The averaged distribution of the intensity (shown 
in both plots by empty circles) has almost no modulation 
with the period of the interference pattern; the observed 
non-uniformity is caused by the interference pattern en- 
velope and the edge effects in the EOM. 

Similar dependencies (with smaller visibilities) were 
measured for the case of pseudo-thermal sources. 



In conclusion, we have demonstrated, both theoreti- 
cally and experimentally, that while the classical visi- 
bility limit for two-photon interference is only 50%, it 
becomes considerably higher in the multi-photon case. 
Three-photon and four-photon interference has been ob- 
served both for coherent and pseudo-thermal light, and 
the maximal visibility values 74% and 93% have been 
achieved. These results, on the one hand, show that clas- 
sical pulsed radiation, in addition to the high spatial res- 
olution of ghost diffraction and ghost imaging, can also 
provide a high visibility, but in this case multi-photon 
correlations should be used. On the other hand, our re- 
sults establish a visibility 'threshold' in multi-photon in- 
terference experiments: for the interference to be clearly 
nonclassical, the visibility should exceed 81.9% in the 
three-photon case and 94.4% in the four-photon one. 
Since exceeding these rather high values is practically 
difhcult, we suggest that the visibility criterion should 
be simply avoided in multi-photon interference experi- 
ments, and some other observable signs of nonclassical- 
ity be used, such as the Lee-Klyshko criterion 18|, re- 
lations between the different-order normalized ICFs 1J| , 
or the scaling of normalized ICFs with the mean pho- 
ton number |19j . This work was supported in part by 
the RFBR grants # 06-02-16393, # 05-02-16391, # 06- 
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FIG. 1: Observation of three-photon interference in the Young 
double-slit scheme. 
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FIG. 2: Experimental setup. GPl and GP2, Glan prisms; S, 
a slit; EOM, electro-optic modulator; BSl, 33%/67% beam- 
splitter; BS2, 50% beamsplitter; Al-3, apertures; PI, 3, glass 
plates; Dl-3, avalanche photodiodes; TCC, triple coincidence 
circuit. 
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FIG. 3: Interference pattern in the normalized third-order 
ICF for (a) coherent sources (b) pseudo-thermal sources ob- 
tained by tilting glass plates at the inputs of detectors 1 and 
3. Empty circles show the intensity distribution given by the 
counting rate Ri of detector 1. Solid lines show the theoretical 
fit with Eq.((2} 




FIG. 4: Single-pulse photograph of the Young interference 
pattern for two coherent sources. The rectangle shows the 
dimensions of the processed area. 




FIG. 5: Interference patterns in normalized third-order (a) 
and fourth-order (b) ICFs for coherent sources obtained by 
processing digital images as shown in Fig. 4. The averaged 
intensity distribution in the x coordinate is shown in both 
plots by empty circles. Thin solid lines show the theoretical 
dependencies ^ (Fig.5a) and Q (Fig.5b). 



